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Abstract: The human arm is kinematically redundant in the task of pointing. As a result, multiple arm 
configurations can be used to complete a pointing task in which the tip of the index finger is brought to a 
preselected point in a 3D space. The authors have developed a four degrees of freedom (DOF)model of the 
human arm with synthesis of two redundancy resolution criteria that were developed as an analytical tool 
for studying the positioning tasks. The two criteria were: (1) minimizing the angular joint displacement 
(Minimal Angular Displacement - MAD) and (2) averaging the limits of the shoulder joint range (Joint 
Range Availability - JRA). As part of the experimental protocol conducted with ten subjects, the kinematics 
of the human arm was acquired with a motion capturing system in a 3D space. The redundant joint angles 
predicted by a equally weighted model synthesizing the MAD and JRA criteria resulted with a linear 
correlation with the experimental data (slope=0.88; ojfset=l°; r 2 =0.52). Given the experiment protocol, 
individual criterion showed weaker correlation with experimental data (MAD slope-0.57, ojfset=14°, 
r=0.36 or JRA slope=0.84, offset=-l°, r=0.45). Solving the inverse kinematics problem of articulated 
redundant serials mechanism such as a human or a robotic arm has applications infields of human-robot 
interaction and wearable robotics, ergonomics, and computer graphics animation. 
Keywords: human arm, redundancy, pointing task, kinematics, optimization. 
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I. Introduction 

Pointing with the fingertip to a preselected point in space is a task that involves three degrees of 
freedom (DOF), which arethe X,Y and Z coordinates of the point in space.The human arm includes seven DOF 
excluding scapular motion. When the wrist joint is fixed, four DOF (^» 7 7'^»^ - Fig. la) remain active. Since 
the number of DOF of the arm is greater than the number of DOF required for the pointing task, the arm is 
considered a redundant manipulator. As such, a specific pointing task can be accomplished by infinite arm 
configurations. As a result, there is not a unique solution for the inverse kinematics (IK) problem involved in 
defining the joint angles of the human arm given a pointing task. 

Despite the human arm redundancy, it has been shown experimentally that a small range of unique 
solutions for the joint angles are selected by human subjects in pointing tasks, a result consistent within and 
across multiple participants [1-5]. It has also been shown experimentally that the final arm configuration depends 
on its initial posture [2-5]. One approach for solving the under-determined IK problem of the redundant human 
arm is by adding additional kinematics, dynamics, or energy-based criteria, formulated as a cost functions. As 
part of the solution, the cost function is either minimized or maximized to provide a unique solution to the IK 
problem when applied to intermediate points along the trajectory of the human arm end effector (i.e. the finger 
tip for a point task) [3- 10]. 

Optimization criteria may be divided into two classes: (1) biomechanical (kinematics and dynamics of 
the human body) criteria and (2) anatomical based criteria. The first class includes the minimal angular 
displacement (MAD) model[3,4], the minimal work model [3,5], the minimal peak kinetic energy model[4], 
the minimal torque change model[5,6] and the minimal potential energy change model [8]. Physical quantities 
such as energy, torque or displacement form the cost function which is further minimized or maximized as part 
of the solution. The second class is based on anatomical models such as the joint range availability (JRA) 
criterion (also called Dexterity) [7] and the minimum discomfort criterion [9]. The cost functions in this class are 
based on anthropometric data of joint motion ranges, the intension of which is to quantify psychophysical 
discomfort related to the proximity to joint limits or nominal arm configuration. 

The majority of these models, when studied individually, and validated experimentally, have 
demonstrated limited capabilities for solving the IK problem of a redundant human arm and predicting arm 
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onfiguration. In order to overcome the limited capabilities of individual criteria, it was suggested that two or 
more criteria should be synthesized with weighted factors [3-5], [12-14]. 



(2) 



II. Materials And Methods 

2.1 Kinematic Arm Model 

The human arm may be modeled as aserial kinematic chain. For the purposes of this study it is modeled 
as a four DOF kinematic linkage, consisting of two links (upper arm and forearm along with the hand) and two 
joints (shoulder joint and elbow joint, wherethe wrist joint is fixed). The shoulder joint is simplified as a ball and 
socket joint with three DOFs, and the elbow is simplified as a revolute joint with a single DOF. This model 
implies that the forearmand hand are aligned with a fixed wrist jointangle during a reach and point movement. 

While the position of the hand in Cartesian space is defined by three coordinates (x,y,z), the posture of 
the arm requires four angles to fully and uniquely specify its configuration. As a result the human arm model 
may be treated as a kinematically redundant mechanism with respect toa pointing task[15]. 

2.1.1 Forward Kinematics 

The forward kinematics equations of a 4 DOF human arm model (see Fig. la and Fig. lb) are defined 
by (1) - (4) for which Sis the pitch angle, Ajis the yaw angle, <fis the upper arm torsion angle (shoulder joint), and 
(fis the flexion/exertion angle of the elbow joint. 

[TH°^]{*?] (1) 
[°fH°^]-ra{7'] 

where °p f =(x f , j f ,z f ) T an d °P e = (jt e , y e , z e ) T are me elbow and fingertip position vectors in a shoulder fixed 
frame, e p f = (ofl-L 2 ) T is the fingertip position vector in an elbow frame, e p f is the elbow's position in a shoulder 
frame, L 2 is the distance from the elbow to the fingertip, °t and s T e are the 4x4 homogenous transformation 
matrices for the shoulder and elbow joints, respectively, as defined by (3) and (4): 

(3) 



(4) 



where T x and T z are 4x4 Homogeneous transformation matrices (pure rotations), T is a 4x4 coordinate 
translation matrix, and Lj is the distance from the shoulder to the elbow. 

Once the position of the hand is fixed at a specific target in space, the elbow joint may swivel around a 
virtual line connecting the shoulder joint and the location of the hand with an angle defined as the swivel angle 
crwhich constitutes the redundancy of the human arm (see Fig. lc). Given the physiological joint limits, the 
elbow joint may only follow a limited arc out of the full circle depicted in Fig. lb. 
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Figure 1:A 4 DOF model of the human arm: (a) Shoulder and elbow joints coordinate systems. (b) The swivel 
angle a. (c) Definition of arm parameters - shoulder joint pitch angle 0, yaw angle 77, torsion angle ( and elbow 

joint flexion/extension angle (p. 

By solving (1) for °p f , we obtain: 

x f = -L^in 77 Sin 0 - Z^Sin ^(Cos^ Sin 77 Cos 0 + Sin £ Cos 77) + Z^Cos <f> Sin 77 Sin 0 (5 ) 

y f = L.CosrjSin 0 + Z^Sin ^(Cos^Sin rjCosO + Sin ^Sin 77)+ L 2 Cos^Sin 77 Sin 0 (6) 
Zf = -L L Co&0 + L z (Sin^Cos^Sin0-Co^Cos<9) (7) 
2.1.2Inverse Kinematics (IK) 

One may note that the FK defines the fingertip position (x,y,z) as a function of the four angles of the arm 6, 
rj, £and 0 For solving the IK, given the fingertip position and the redundancy of the mechanism, the problem is 
underdetermined. There are 4 unknowns (0, A], <f and (p) and only 3 equations (5-7). To obtain a closed form 
solution one of the unknowns must be specified. Alternatively, the IK problem can be formulated by using the 
swivel angle (8) or by defining the elbow joint position (9). These two alternative approaches were used to 
solve the IK namely the arm configuration 0=(0,/7,£<$ T : 

0 = (P f ,a) (8) 

& = (Pf,P.) (9) 
wherep^^^zj 



F 

(fingertip) 




Figure 2: Definition of a local coordinate system at the center of the swivel circle, the point Q and the angle y. 

In order to find the arm's posture, given a swivel angle, first, the elbow's position is defined first. A 
virtual line is defined connecting the center of the shoulder joint and the fingertip. This line is depicted in Fig. 2 
as the swivel axis OF A plane perpendicular to this line intersects with the line a point C and with the center of 

the elbow joint at point E. The swivel circle is formed at this plane with a center at point C with a radius r. A 
local coordinate system is define with the origin at the center of the circle at point C and three unit vectors u,v, 
and n such that uis pointing down along the vector of gravity, v is located at the same plane but perpendicular to 
u, and n is perpendicular to the plane. The three unit vectors are defined in mathematically as : 



11= - 



(10) 



U = 1 



-z + ( z • n )n 



(11) 



||-z + (z-n>i|| 

v=uxn (12) 
where z is a unit vector with the direction of the z axis as defined in Fig. la. We also define the angle y = ZFOE 
and calculate it as in: 



y = Cos 



||of|| 



4 



2L t L 2 



The center of the swivel circle is defined by: 
OC = nl^Cosy 



(13) 



(14) 
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and the radius of the circle is defined by: 



Z?- OC 

The elbow's position P e = OEcan now be obtained by: 
OE = OC + r(uCosa + vSina) 

The 4 DOF angles are calculated using trigonometric relations as follows: 



e = cos 1 



77 = Atan2(-x e ,.y e ) 

G = Atan2(L 1 (x e y f - x t y e ), y e (y e z, -y f z e )-x e (z e x f - z,x e )) 



> = 7t - COS 



To find the arm's posture, given an elbow position, one may notice that: 

OQ=OC+ra 

and then: 

CE v ^ ! CE CQ 
a = —= — ^Cos 



CE-v 



CE-CQ 



(15) 

(16) 

(17) 

(18) 
(19) 

(20) 

(21) 
(22) 



Note that none of these procedures provide a solution to the IK problem of the kinematically redundant 
mechanisms, because either swivel angle, or elbow position is required. 



2.2 Shoulder and elbow joints motion range 

An experimental model was previously developed which predicts the shoulder complex range of 
motion [16]. This model defines the shoulder sinus cone, which restricts the angular motion of the shoulder 
joint's pitch and yaw angles as depicted in Fig. 3a. 

The humeral torsionmotion range was found to be dependent on the former two angles depicted in Fig 
3b using a model based on experimental data [17]. The maximal internal and external upper arm torsion surfaces 
were fitted into a polynomial function. Fig. 4 shows the upper and lower limits of the torsion range of 
motion.The motion range of the elbow joint (flexion/extension angle) is bounded by a minimal and a maximal 
value as defined by (j) m ^ n <(/)< ^ max . Since the elbow flexion/extension angle is uniquely defined by the distance 
between the center of the shoulder joint and the hand, it can be calculated directly according to the hand's 
position and the segments' lengths. 

(a) (b) 

200 



V 




-150 



Shoulder yaw // (cleg) 




Figure 3: (a) Shoulder joints' dependencies on pitch and yaw angles, (b) Upper (red) and lower (green) limits of 
humeral torsion. The distance between the two surfaces defines the torsion motion range. 

Using this model, we can define a valid posture of the human arm to be a posture that complies with 
the conditions in (23): 



^ mi n(^)<^<^ m ax(^7) 



(23) 
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An alternative and simpler approach to the shoulder's joint motion range is given by constant ranges to each 
DOF, independent on each other. Values for joint limits can be found in anthropometric data sources [18] as in 
(24): 



(24) 



f0<6><180° 
-134° < rj < 48° 
-34° < £ < 97° 

Both approaches are used in this study, in order to compare their effect on the solution of the kinematic 
redundancy problem of the human arm. 



2.3 Joint Range Availability (JRA) Criterion 

The JRA criterion is based on the assumption that the human arm tends to adopt postures with joint 
angles that are as close as possible to their mid-range values and as far as possible from their joint limits. As the 
elbow joint swivels, and (theoretically) provides an infinite number of possible arm postures (solutions to the IK 
problem), the pitch, yaw and torsion angles of the shoulder joint are adjusted appropriately to maintain the hand 
position. A valid anatomical posture is achieved if all three angles of the shoulder joint are within their 
anatomical ranges of motion either according to (23) or according to (24).There is a continuous subset of valid 
arm postures with a swivel angle in the range of a <a<a • Based on the JRA criterion, the optimal posture 

mm max 

is achieved by a mean value of the swivel angle limits defined by: 

(25) 



mm max 



Figure 4 illustrates the swivel angle limits and the optimal posture at some handposition in space. 

A Z 




ingertip 



elbow 



Figure 4: Swivel angle limits and optimal posture. The solid arc denotes valid postures domain. Dashed arc 

denotes invalid postures domain. 
2.4Minimal Angular Displacement Criterion (MAD) 

The minimal angular displacement (MAD) criterion minimizes the sum of the differences of the 
various joint angles, between their initial and final values. In other words, the final arm posture defined by this 
criterion yields the shortest distance between the initial and the final value of the joint angles in joint space. This 
criterion can be formulated as an optimization problem using the following cost function: 

/mad = (0* ~ # f + [rf -jf +(<?-<? f (26) 

where i denotes initial value and f denotes final value for the pitch, yaw, and torsion angles. However, it has 
been shown experimentally that the shoulder joint velocities vary simultaneously [4]. Therefore, the shoulder 
shifts directly from its initial pitch and yaw angles to its final pitch and yaw angles (see Fig. 5), and not by two 
separate rotations as implied by (26). This rotation occurs about an axis that is normal to the plane formed by the 
shoulder (O), the initial elbow position (Pi), and the elbow final position (P 2 ). The direct rotation of the 
shoulder, and the direct rotation angle are illustrated in the following figure: 
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Figure 5: The direct rotation angle 8. Pi and P 2 are two elbow positions, and m is a vector normal to the plane 

defined by points O, Pi and P 2 . 

The direct rotation angle is calculated by: 



^Cos 1 



P P 

1 1 r 2 

Hp II -Hp II 

1 2 



and: 



Therefore, the cost function of MAD criterion in this study is formulated as follows: 

\2 



f, 



MAD 



=s 2 + 



(27) 



(28) 



(29) 



Since a 4 DOF arm model is redundant by only one DOF, this optimization problem can be solved by a one 
degree step grid search where solving for the swivel angle that minimizes the cost function under the given 
constraint, using (29). During the brute force search, the initial posture of the arm remains constant, while the 
final posture varies with the value of the swivel angle, and the cost function's value changes accordingly. 



2.5 Bi-Criterion Model 

The JRA and MAD criteria are used independently, as previously explained, to calculate the swivel 
angles CT JRA and Cr MA D,respectively. The synthesized criterion model then merges the two results by calculating 
the weighted average of the swivel angle OC optiYml as defined by: 

^optimal = ^MAD + ( l " ^KrA (30) 

where kis a weight factor. Its value is later optimized to match experimental results with the model prediction, 
and to explore its effect on the correlation with the experimental results. 

2.6 Pointing Experiment 

2.6.1 Subjects Definition 

Ten right handed volunteer (ages 22 years to 28 years, heights 1.65 m to 1.93 m, weights 52 Kg to 94 
Kg) performed unconstrained point to point movements. None of the participants reported having any clinical 
symptoms or any history of motor, sensory, or neurological disorders. The subjects were naive as to the purpose 
the experiment, and were only instructed to perform point to point reaching movement in a natural self -paced 
way. The experimental protocol was approved by the UCSC IRB committee. 

2.6.2 Experimental Setup 

Subjects sat on a chair with back rest and adjustable height. They were positioned in front of a 
spherical target array (see Fig. 6a and Fig. 6b) with a radius of 750 mm, such that the right shoulder was in line 
with the center of the sphere.The distance between the chair and the target was adjusted for each subject so that 
his/her elbow flexion angle was approximately 90°. 
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Figure 6: (a) The spherical target array. The stars denote 44 available targets, (b) A subject pointing to a target. 

Out of the 44 available targets, eight targets were selected (Fig. 7Figure 7), out of which four (1, 4, 5, 
and 8) are at boundary of arm's reaching workspace, whereas the remaining four targets are closer to the center 
of the workspace (2, 3, 6, and 7). Targets consisted of bar ends, retained by a frame, specifically designed for 
pointing experiments. 




Figure 7: The eight targets used in the pointing experiment, denoted 1 through 8 (large figures). 



2.6.3 Experiment Protocol and Data Collection 

As part of the experimental protocol the subjects conducted eight sessions of pointing movements with 
their right arm. Each of the eight sessions consisted of five sets of movements. In each set, subjects started 
from each one of the seven targets and moved their fingertip to the selected destination target of that set, thus 
performing seven movements. Fatigue was avoided by resting period between each session. Each subject 
conducted a total of 280 pointing tasks. 

Reflected markers were attached to the human arm and upper body using a standard model. Position of 
the markers were acquired by a 10 cameras motion capture system (Vicon Motion Systems) at a sampling rate 
of 100Hz. From experimental sessions, the initial and the final postures of the human arm were identified. Joint 
angles were calculated off line according to (17)-(20). For each subject, 56 average initial and final postures 
were calculated. 



III. Results 



3.1 Experimental Results 

Postural invariance of human arm pointing movement was demonstrated in the experiment. Each 
pointing movement was repeated 5 times. The swivel angle associated with repeated final pointing postures 
varied in an average range of 8.1° with a standard deviation of 4.5°. This variance is relatively small compared 
to the variance of average final posture when pointing to all eight targets (68° for all subjects). Fig. 8 illustrates 
the variation of final posture and the effect of initial posture on the final posture for subject 1. 
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Task Number 

Figure 8: Final posture variability for subject 1. 



3.2 MAD model 

The correlation between the experimental results and MAD prediction model, using equations (26)is 
depicted in Fig. 9a. Its trendline slope is 0.5, the offset is 16°, and r 2 =0.37. The correlation between the 
experimental results and MAD prediction model, using equation (29) is depicted in Fig. 9b. Its trendline slope is 
0.57, the offset is 14°, and r 2 =0.36. These two versions of the MAD model achieved similar results with only a 
slight advantage for the case of using equation (29). 




predicted swivel [deg] predicted swivel [deg] 



Figure 9: Correlation between experimental results and predictions of MAD model using (26) (a) and (29) (b) 

for all subjects. 

3.3 JRA model 

Two alternative definitions for the joint limits were previously defined by equations (23) and (24). 
Using each of these joint limits resulted in two different arm configurations and associated with two different 
swivel angles. A typical result expressing these differences is depicted in Fig. 10. The average difference 
between the two joint limits is 26.2° with a standard deviation of 14.4°. 
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Figure 10: Comparing swivel angle limits according to two different approaches, wheredots represent the 
targets,thick lines representupper arms,thin lines represent fore-arms,blue lines represent the extremal valid 
postures found according to (23),green lines represent the extremal valid postures found according to (24),the 
arc represents the trajectory of the elbow between swivel extremes, andthe numbers in the brackets indicate the 
differences between the two approaches for pointing to target 2. 

The variation in final swivel angle for the spherical surface passing through the targets is depicted in 
Fig. 11. Uncolored areas in Fig. llrepresent unreachable positions due to either being too far from the shoulder 
(compared to the total arm length), or requiring an elbow flexion angle greater than 150° which occurs in the 
vicinity of the shoulder. Occlusion of the fingertip with the body is ignored in this simulation. 




Figure 1 1 : Variation in final swivel angle predicted by the JRA criterion as a function of target position for the 

spherical surface containing the targets (black dots). 

The correlation between the experimental results and JRA prediction model, using equations (23) has a 
trendline slope of 0.91, an offset of -28° and r 2 =0.27, whereas using equation (24) showed a stronger correlation, 
with a trendline slope of 0.84, an offset of -1° and r 2 =0.45. These correlations are depicted in Fig. 12a andFig. 
12b, Respectively. 




predicted swivel [deg] 



10 20 30 40 50 60 70 

predicted swivel [deg] 



Figure 12: Correlation between experimental results and Predictions of JRA model for all subjects (a) using 

(23), and (b) using (24). 
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3.4 Bi-criterion model 

The bi-criterion model requires a selection for the weight factor k. Three approaches for selecting the 
weight factor k: (1) objective sum, in which k=0.5, as suggested in[20][20]; (2) optimizing k for highest value 
of r 2 ; and (3) optimizing k for best trendline slope of unity. 

For the bi-criterion model, MAD and the JRA criteria utilized equations (29) and (24) respectively. The 
correlation between the experimental results and bi-criterion model predictions with a weight factor k=0.5, has a 
trendline slope of 0.88, an offset of 1° and r 2 =0.52 as depicted in Fig. 13. 



£ 50 



CD 40 



CD 

E 



slope = 0.88 
offset = 1 5 
r 2 = 0.52 






sub 1 




sub 2 




sub 3 


+ 


sub 4 


+ 


sub 5 


□ 


sub 6 




sub 7 




sub 8 




sub 9 


< 


sub 10 





10 20 30 40 50 60 70 

predicted swivel [deg] 



Figure 13: Correlation between experimental results and Predictions of the combined model, where k=0.5 

(objective sum), for all subjects. 

The effect of the weight factor value (k) on the trendline slope and on r 2 is depicted in Fig. 14. The 
highest value of slope (0.29) was achieved with &=0.29, and the highest value of r 2 =0.54 was achieved with 
k=035. Utilizing the objective sum approach (k=0.5) yielded a slope of 0.88 and r 2 =0.82. 
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Figure 14: The effect of the weight factor value (k) on the trendline slope (dashdot) and on r (solid line) for the 
bi-criterion model combining MAD and JRA criteria. The numbers in parentheses are k and slope or r 2 . 

IV. Discussion 

The focal point of this reported research effort is modeling synthesis approach for redundancy 
resolution of the human arm in pointing tasks. Results for the single criterion models (MAD and JRA) and for 
the bi-criterion model are presented in table I: 
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Table I: Single and bi-criterion models results 



Model 


Method 


K 


Slope 


offset 


r 2 


MAD 


Eq. (26) 




0.5 


16° 


0.37 


MAD 


Eq. (29) 




0.57 


14° 


0.36 


JRA 


Eq. (24) 




0.84 


-1° 


0.45 


JRA 


Eq. (23) 




0.91 


-28° 


0.27 


MAD+JRA 


Objective sum 


0.5 


0.88 


1° 


0.52 


MAD+JRA 


Optimal slope 


0.29 


0.94 


-3° 


0.53 


MAD+JRA 


Optimal r 2 


0.35 


0.92 


-2° 


0.54 



For a 4 DOF model of the human arm, where there is only one redundant DOF, the synthesized model 
using two criteria provided higher correlation with our experimental data, compared to results achieved by two 
single criterion models. 

Utilizing the MAD criterion with the direct rotation angle using (29) showed only a minor 
improvement in the correlation, compared to utilizing it with (26). Swivel angle limits that were calculated 
according to (23) were found to differ significantly (average 26.2°) from those calculated according to (24). 
However, using the more realistic description of the shoulder motion range [18], where the joint limits are co- 
dependent, did not improve the correlation given the current experimental data. Since the correlation of the JRA 
model using (23) had a large trendline offset (-28°), this result is not surprising. The large trendline offset may 
be associated with fact that (23) allows a larger medial torsion range (up to 97°) than (24) (up to 220°), as 
illustrated in Fig. 15. 

50 

34 




-100 -50 
n [deg] 

Figure 15: Comparison between humeral torsion range according to (23) (surfaces) and according to (24) 

(straight lines). 



Although the JRA model, using (24) had a relatively close to unity trendline slope, and a relatively high 
r 2 , it cannot be considered to be a good prediction model due to the fact that it is a posture based model, which 
does not account for the effect of initial posture on the final posture. Other single criterion models shown in 
Table I had relatively low slope or low r 2 . Therefore, it may be deduced that synthesizing the two criteria 
achieves better correlation given the current experimental data, compared to the correlation with the predictions 
of single criteria (MAD or JRA) even with the objective sum approach, which does not optimize the weight 
factor to match the experimental data. 
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Resolution of human arm redundancy in point tasks by synthesizing two criteria 



Previous research efforts indicated that correlation of alternative criteria such as the minimal work and 
minimal torque change were associated with a slope of 0.3 and r 2 =0.56 for both criteria [5]. Using the minimum 
peak kinetic energy criterion[4] led to correlation with slope of approximately 1 and r 2 in the range of 0.522 to 
0.915 using a database of four subjects. 

Future work will focus on examining the various approaches for integrating criteria (objective sum, 
min-max, and global sum), and on finding an optimal combination of several criteria to achieve good correlation 
with experimental data. 

Redundancy resolution of articulated serial mechanism such as human or robotic arm has applications 
in the field of wearable robotics [21,22], where the robot and human body are physically coupled and the need 
to implement a redundancy resolution algorithm into the control system of the robot is essential to guarantee 
natural integration between the two systems. Additional application that may benefit for a solution the problem 
under study might be simulations in ergonomics and computer graphics. 
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